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We investigate thermoelectric transport through Andreev interferometers. We show that the ratio 
of the thermal and the charge conductance exhibits large oscillations with the phase difference <j> 
between the two superconducting contacts, and that the Wiedemann-Franz law holds only when 
(f> = 7T. A large average thermopower furthermore emerges whenever there is an asymmetry in the 
dwell times to reach the superconducting contacts. When this is the case, the thermopower is odd 
in <j>. In contrast, when the average times to reach either superconducting contact are the same, 
the average thermopower is zero, however mesoscopic effects (analogous to universal conductance 
fluctuations) lead to a sample-dependent thermopower which is systematically even in (j>. 

PACS numbers: 74.45.+C, 74.78.Na, 73.23.-b 



Introduction. The processing of information un- 
avoidably generates heat [l|. Conventional micro- 
electronics use electric potentials to switch currents on 
and off. At the nanoscale, however, this becomes energet- 
ically prohibitive and generates an amount of heat that is 
hard to dissipate. As new architectures are explored for 
quantum communication and computing, the question of 
dissipating heat is again of central importance. In cur- 
rent prototypes for quantum information processors, one 
of the slowest steps is cooling down the qubits in be- 
tween computations Q, thus it is crucial to understand 
heat flows at the nanoscale and sub-Kelvin temperatures 
where quantum coherent effects are ubiquitous. 

While quantum interference effects occur in all meso- 
scopic systems |3( , many of them are hugely magnified by 
Andreev reflection in normal- metallic/superconducting 
nanostructures [1, H, 0, 0] • Experiments on Andreev in- 
terferometers - metallic constrictions contacted to two 
superconducting terminals with a phase difference, </> - 
have revealed thermoelectric properties that are strongly 
affecte d by these magnified interference effects [jj H, MM, 
El, El El El- The device properties can be probed (and 
controlled) by varying <fi with an applied magnetic flux 
or a supercurrent. For instance, the charge, G, and ther- 
mal, 2, conductances and the thermopower, S, oscillate 
periodically with </>. The salient observations are that (i) 
the amplitude of the conductance oscillations can largely 
exceeds e 2 jh and is typically larger in samples with larger 
average conductance, (ii) G has its maxima where S has 
its minima and vice- versa in violation of the Wicdcmann- 
Franz law, (iii) S is significantly larger than in normal 
metals in absence of superconductivity, (iv) S is either 
even or odd in <f>, depending on the interferometer geom- 
etry, and (v) S exhibits oscillations of maximal amplitude 
at an intermediate temperature. Despite extensive the- 
oretical investigations 0, El, El, E3, El, El, El, EI , a 
unified theoretical picture of these observations is still 
lacking. In particular, the existing scenarios for an odd 
thermopower are always associated with a temperature 



gradient between the two contacts to the superconduct- 
ing terminals 0, 0, El. it is thus unclear whether 
these theories can capture the recently observed odd ther- 
mopower [ll| in the geometry of Fig. [TJl, where this 
gradient most likely vanishes. 

Motivated by these experimental findings, we inves- 
tigate experimentally relevant models of Andreev inter- 
ferometers where ideal metallic leads carrying Ni 3> 1 
modes are connected at either chaotic ballistic or disor- 
dered quantum dots with no spatial symmetry. The dots 
are contacted to two s-wave superconductors with order 
parameters Ae 1 ^ , each carrying iVg j channels, i = L, R. 
Physical properties depend only on the phase difference 
4>h — 0r, so we set (f>h = 0/2 and 0r = — </>/2 with A 6 H. 
We take the superconductors to be islands through which 
no current flows on time average in steady-state, as at 
propriate to the experiments of Rcfs. 
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The models are sketched in Fig. Q] and are devised to 
have the same topology as the house (a and b) and par- 



allelogram (c) interferometers of Refs. [liJ, [ll||, and (d) 
the interferometer of Ref. (l3j which we call hot-middle. 
These models differ by the absence (a) or presence (b,c,d) 
of correlation between the action phase a quasiparticlc 
accumulates and the superconducting phase it acquires 
at Andreev reflections. This correlation is key to under- 
standing the large thermopowers induced by the presence 
of superconductivity, because it breaks particle-hole sym- 
metry. 

The theory we are about to present gives a unified pic- 
ture of thermoelectric transport through mesoscopic An- 
dreev interferometers. Extrapolated to diffusive systems, 
it exhibits all the main experimental observations listed 
above. For the asymmetric house (model b), parallelo- 
gram and hot-middle interferometers, it predicts that on 
average the thermopower will be an odd oscillatory func- 
tion of <$>, even in the absence of a temperature difference 
between the two superconducting contacts. In contrast, 
for the symmetric-house interferometer it predicts that 
the average thermopower is zero, but that mesoscopic 
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(a) symmetric house (b) asymmetric house (c) parallelogram (d) hot-centred C 




Figure 1: The four Andreev interferometers considered in this letter, (a)-(b) Single-cavity, two-terminal models, topologically 
equivalent to the house interferometer of Ref. [13] ; the symmetry between the leads to the superconductin g co ntacts is broken 
in model (b). (c) Double-dot model, topologically equivalent to the parallelogram interferometer of Ref. [1Q |. (d) Triple-dot 
model topologically equivalent to the hot-middle interferometer of Ref. []]§], with a hotter lead M. 



fluctuations render the thermopower random in sign, but 
systematically even in <f>. 

Thermoelectric transport. The linear response ex- 
pression for charge, Ii, and heat, J^, currents in lead i (l7j 
(summation over doubly-occurring indices is assumed) is 



with the derivative F'(e) of the Fermi function and the 
base temperature T. Since there is no net current into the 
S loop, its potential, Vq, is tuned to ensure J2j Ij = 0- I n 
two-terminal geometries, Fig.QJi-c, Eq. ((T|) reduces to I = 
G(V h -V K )+B(T h -T R ) and J = r(V L -F R )+S(T L -T R ). 

Taking T"- (e) as the transmission coefficient for a j3- 
quasiparticle (e or h) injected from lead j at energy e 



exiting as an a-quasiparticle in lead i, Ref. 17J gives 



Gtj (e) = \ [2NAJ -T^ + T^ + T^-T^], (2a) 
Eij (e) = ^ -T%-Tg-1tt- T^] > ( 2b ) 



Iff 

Bij(s)= -—[T^-T^ + T^-T* 11 ], 
Ipf 

r«(e) = -—[r^ + Tt-Ti 
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(2d) 



Average diagonal thermoelectric coefficients. 

To evaluate the transmission probabilities for Eqs. (J5J), 
we use the Feynman rules in Ref. 22j, taking into ac- 



count the presence of two S contacts with phase differ- 
ence (j>. We work perturbativcly in the ratio Ns/N of 
the total number of channels carried by the S contacts 
to those carried by the normal leads. All contributions 
up to Q[(Ns/N) 2 } are shown in Fig. [21 the </>-dependent 
ones are ee2II, he2I, he2II and he2III. For model a and 
b, their value in Eqs. (3) and (4) of Ref. [22| is now mul- 
tiplied by 5 (0) - [Af L + Af R + 2N sh N SR cos 0]/JV|. For 
Nl = N R = N, the average charge and thermal conduc- 
tances (neglecting weak-localization) are 



(G) = N/2 + Nlg(<f>)f G (T)/8N, (3a) 
(B> = f^L(N/2-Nlg(cb)f E (T)/8N). (3b) 



The main difference between charge and heat con- 
ductances in this symmetric configuration is that the 
backscattering contribution he21 in is absent in Tp£, 
which brings the periodic oscillations in G and S out of 
phase by 7r, with G bein g m inimal at <f> = 0. This fits with 
the experiment of Ref. lf| . The thermal dampings are 



polynomial, generalized zeta and poly gamma functions 
of a = AkBTTn, where td is the dwell time. Asymptoti- 
cally / G (0) = 1, f 3 (a < 1) « 10.4a, and f G (a) w n/a, 
/a (a) ~ 32/7r 2 a, for a > 1. 

Eqs. imply that there are coherent oscillations of 
the Wiedemann-Franz (WF) ratio 
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Nlg(<f>)[f G (T) + f s (T)] 



4N< 



,(4) 



where our results below show that TS is small enough to 
neglect, and Iq = 7r 2 fc|/3e 2 is the Lorenz number. Thus, 
unless 4> = 7r and ATsl = Na, R simultaneously, supercon- 
ductivity causes a 0[(Ns/N) 2 ] violation of the WF law, 
parametrically larger than the OfA^ 1 ] violation induced 
by mesoscopic fluctuations in metallic samples (23|. 

Average thermopower. The off-diagonal thermo- 
electric coefficients satisfy the Onsager relation B = 
— r/T. We checked that our theory preserves this sym- 
metry and only discuss B from now on. In two-terminal 
arrangements, G is symmetric in <j). Therefore the sym- 
metry of the thermopower coefficient S = —B/G is deter- 
mined by the symmetry of B. To leading order in A^, par- 
ticle hole symmetry, e —> —s, is equivalent to reversing 
the superconducting phases, (f> — > — 4>. Combining this 
with Eq. (f2"c|) , we straightforwardly conclude that (B) is 
generically odd in 0, up to weak-localization corrections. 
In a house geometry, however, one can interchange the 
superconducting leads, and thus reverse the supercon- 
ducting phases without changing the physics. Thus (B) a 
must be even in <fi. Neglecting weak localization correc- 
tions, one thus has (B) = and S = for model a. 

How can a finite leading-order thermopower emerge? 
Our symmetry argument breaks down when there arc 
correlations between the action phase a quasiparticle ac- 
cumulates on its way through the system and the su- 
perconducting phase that it picks at Andreev reflections. 
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Figure 2: (Color online) Contributions to (T^j) (first three) and (T^f) (last four). Green (violet) paths indicate electrons 
(holes), dashed lines indicate complex-conjugated amplitudes and circles show trajectory encounters. Normal leads are labelled 
z,j while superconductors are labelled Sa, S/3. Contributions to (T* h ) and (T?f) have e and h interchanged everywhere. 



We find that, when present, these correlations generate 
a finite, odd average thermopower. This is most eas- 
ily seen by analyzing the asymmetric house interferom- 
eter of Fig. [TJd, where a neck renders the journey to- 
ward SR systematically longer. For simplicity, we as- 
sume that all trajectories going into the neck spend a 
time St in it before they hit the SR lead. Now take con- 
tribution he2I. If the solid path hits SR and the dashed 
hits SL, it induces a phase of <p from the S leads, and 
a phase of 2eSr from the extra length of the solid path. 
If the solid path goes to SL and the dashed to SR, we 
get the opposite phases. Interchanging e and h, means 
4> — ► — </>, thus these contributions to T^ e — Tf- 1 behave 
like cos(2e^r + 0) — cos(2e<5r — 0). The prefactor on this 
contribution is easily found using the Feynman rules in 
Ref. [2^. Treating the other contributions in the same 
way, we find that the leading-order average thermopower 
for model b is 



(B h ) 4fc B iVsLiVgR 



<G b ) e (N L + N R y 



I b (T) sin< 



(5) 



where I b (T) = -(/crT)- 1 dee F'(e) sin(2e<5r)/(l + 
4e 2 t d ). This energy integral is zero for St — (sym- 
metric house), is linear in T for T <C St" 1 ,t i J , is max- 
imal when fceT ~ St" 1 ~ r^ 1 and decays as T" 1 when 
T > t d 1 for r D > St, and as exp[-27rfc B TciT] for T > St 
for St 3> td. The average thermopower is always odd in 
(j>, but we stress that for it to be finite we need a system- 
atic asymmetry in the distributions of path lengths to SL 
and SR. An asymmetry in the probability of hitting the 
two S-contacts, such as for A^l 7^ Nsr, is not sufficient. 

The presence of the neck with N n channels, in the par- 
allelogram interferometer of Fig.[T}; also breaks symmetry 
between the length of paths to SL and paths to SR. Wc 
consider N D <C N, and treat the problem to leading order 
in N n /N. The two cavities are not symmetric, and might 
have different dwell times tdl and tdr. To leading order 
in N n /N and N s /N we obtain 



(5 C 



2fc B N SL N SR N n (N SL - N SR ) I C (T) sin^ 
e N L N R [N n (N SL + N SR ) + 2N SL N< 



SB 



(6) 



where J C (T) = -(knT)- 1 J^°de e F'(e)2Im[c i2eST /A(e)] 
with A(e) = (1 + i2£Tdl)(1 + i2er DR ). Again the ther- 
mopower is odd in but to be finite it requires A^sl 7^ 



iVsR. In fact, for iVgL = Ng R there are odd-0 con- 
tributions to (S c ) at next order in Ns, n /N (assuming 
tdl 7^ tdr). We note in passing that a similar expres- 
sion is obtained for the hook geometry of Ref. [l(j ■ 

We apply this (7V n /iV)-perturbation theory to the 
three-dot model of Fig. [T}i. The L and R leads have 
voltages Vl,Vr such that no current flows in any lead 
(or S contact) when lead M is held at a temperature 
T + ST. The thermopower S% M = V a /ST, a <E L, R, is 
lead-dependent. To find this, we solve Eq. (|TJ) for Vb, Vl 
and V R when all charge currents are zero. One obtains 



<S d LM > = 



2k T 



N? lR N nh N Sh N SR I d (T) sine/. 



e N M N L N R [N nL N nR + N' s (N nL + N nR )] 



,(7) 



where N' s 



2N SL N SR /(N sh + N SR ), and J d (T) is given 
by I C {T) with St — ► (Stt_, + St r ) and A(e) gaining a 
factor of (1 + 2ieT DM )- Similarly we find (Sf M ) = 
-(N nh /N nR )(S% M ). Thus for N nR ~ A nL , one obtains 
(Sf M ) ~ -(S^ M ). All this fits with the experimental 
data or Ref. [13j . 

It is worth recalling that Eqs. ([5HZ]) are leading order 
in Ns/N, and thus neglect oscillations cx Ng cos (f>/N in 
the denominator of (S) (coming from (G)). These terms 
generate higher odd-tfi harmonics, not unlike the experi- 
mental findings. 

Mesoscopic fluctuations. Since the average ther- 
mopower vanishes for model a, we look at mesoscopic 
fluctuations. We consider contributions to (Tp£ — Trl + 
7r£ — T^'l) before mesoscopic average. The diagrams 
shown in Fig. [3] give a contribution (Tp£ — T$£) oc 
sin[e(T-t 2 )] sin[i2E F (T+t 2 ) + (p} for T = ti+i 3 — *4- Wc 
sum over the 24 permutations of the four trajectory du- 
rations ti, i = 1, 2, 3, 4, which generally gives both even-<fi 
and odd-0 contributions. For the special case of the sym- 
metric house, however, the odd contributions cancel out 
exactly, since for every contribution touching both super- 
conducting contacts, there is a contribution with equal 
weight touching the superconducting contacts in the re- 
versed sequence. These two come with opposite signs of 
0, so the sum is even in <f>. An analysis of Xr£— Trl shows 
the same behavior, therefore the sample-dependent ther- 
mopower must be even in <f> in this case. 

To evaluate the typical magnitude St yp of the ther- 
mopower in a single measurement, we estimate for var_B a . 
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Figure 3: Contributions to ee2II and hh2II before ensemble 
averaging. 



fluctuations. 
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This is done by pairing any two contributions in Fig. [21 
which necessitates to add at least two encounters. From 
the Feynman rules in Ref. (22|, one obtains the leading 
order in Nq/Nt by pairing eeO with eeO, ee2II and he2III, 
including all permutations with e<-> h. Then B typ ~ 
rms[B a ] ~ (ir 2 klT/3e)T B [ Kl + K 2 N SL N SR cos((b)/N 2 ] at 
low temperature AibTYd <C 1 where Ki t 2 = 0[1]. The 
first term is dominated by normal metal fluctuations of 
B 24| . In contrast to Aharonov-Bohm oscillations, meso- 
scopic fluctuations of thermoelectric coefficients are not 
paramctrically enhanced by the presence of superconduc- 
tivity. A single measurement of the house interferome- 
ter of Ref. [10| thus typically produces an even— ther- 
mopower, 



Si 



typ 



(2k B /e)(S + S4A SL A SR cos(0)/A 2 



(8) 



which need not vanish at <f> = 0. The constants So.i ~ 
(2e 2 //i)G _1 <C 1, so the above cos^-term is in princi- 
ple smaller than the odd thermopowers, Eqs. (JHHl]). For 
a given sample, 6*0,1 are random in sign. In a sym- 
metric house geometry, Ref. [Iol | reported an even ther- 
mopower at T = 38mK of similar magnitude as the 
odd thermopowcr found in a parallelogram geometry at 
T = 350mK. This might be due to a much stronger tem- 
perature damping in the latter case, or to Asl ~ Asr in 
the parallelogram, or both. We note that the experiments 
found S a (cf) = 0) 7^ 0, in agreement with our theory. 

Conclusions. Our theory potentially explains all ex- 
isting thermopower experiments on Andreev interferom- 
eters 0, H, [lfl El B El • Finite average thermopowers 
are systematically odd in </>, and emerge when the geom- 
etry correlates trajectory durations to superconducting 
contacts with the phase at the contacts. When the aver- 
age thermopower vanishes, mcsoscopic fluctuations can- 
not be neglected. We discussed the latter for the first 
time and showed that t hey are systematically even in <fi. 
Unlike earlier theories 0, 0, EH, our mechanisms for 
thermopower do not presuppose charge imbalance nor re- 
quire temperature differences between the superconduct- 
ing contacts. Waiving the latter requirement allows us 
in particular to explain the odd thermopower recently 
found in hot-middle interferometers [HI]. We hope that 
the validity of our theory for even S in house interferom- 
eters will soon be checked by investigations of mesoscopic 
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